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Abstract 

A theory of Electric Dipole Spin Resonance (EDSR), that is caused by various mechanisms of 
spin-orbit coupling, is developed as applied to free electrons in a parabolic quantum well. Choosing 
a parabolic shape of the well has allowed us to find explicit expressions for the EDSR intensity and 
its dependence on the magnetic field direction in terms of the basic parameters of the Hamiltonian. 
By using these expressions, we have investigated and compared the effect of specific mechanisms of 
spin orbit (SO) coupling and different polarizations of ac electric field on the intensity of EDSR. It 
is our basic assumption that the SO coupling energy is small compared with all different compet- 
ing energies (the confinement energy, and the cyclotron and Zeeman energies) that allowed us to 
describe all SO coupling mechanisms in the framework of the same general approach. For this pur- 
pose, we have developed an operator formalism for calculating matrix elements of the transitions 
between different quantum levels. To make these calculations efficient enough and to derive explicit 
and concise expressions for the EDSR intensity, we have established a set of remarkable identities 
relating the eigenfrequencies and the angles defining the spatial orientation of the quantizing mag- 
netic field B (6, 4>) . Applicability of these identities is not restricted by EDSR and we expect 
them to be useful for the general theory of parabolic quantum wells. The angular dependences of 
the EDSR intensity, found for various SO coupling mechanisms, show a fine structure consisting of 
alternating up- and down-cusps originating from repopulating different quantum levels and their 
spin sublevels. Angular dependences of the EDSR intensity are indicative of the relative contribu- 
tions of the competing mechanisms of SO coupling. Our results show that electrical manipulating 
electron spins in quantum wells is generally highly efficient, especially by an in-plane ac electric 
field. 

PACS numbers: 71.70.Ej, 76.20. +q, 78.67.De, 85.75.-d 
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I. INTRODUCTION 



Efficient manipulation of electron spins by an external ac field is one of the central 
problems of semiconductor spintronics,— & quantum computing 3 and information processing.— 
The original proposals regarding spin manipulation were based on using a time dependent 
magnetic field H(t) . However, there is a growing understanding of the advantages of spin 
manipulation by a time- dependent electric field E(t) that allows the access to electron spins 
with nanometer precision and can provide much stronger coupling to electron spins through 
various mechanisms of SO interaction.— £ Different options of electrical manipulating electron 
spins in semiconductor nanostructures range from adiabatic pumping spin currents from 
quantum dots^ to manipulating electron spins in quantum wells (QW) at the spin resonance 
frequency by employing various mechanisms of SO coupling. Spin orbit interaction 

Hso = H OTh (k, <r) + H z (r, cr) (1) 

can be usually represented as a sum of the orbital contribution H OT b( k , cr ) depending on the 
momenta k and Pauli matrices cr , and the Zeeman contribution Hz(r , cr) depending 
on the coordinates r and matrices cr . Because both terms include orbital operators ( k 
or r ) and matrices cr , they represent different mechanisms of SO coupling. 

For two-dimensional (2D) electrons in QWs, two basic mechanisms of the orbital SO 
coupling are directly related to the QW symmetry properties. They stem from the structure 
inversion asymmetry (SIA) mechanism described by the Rashba term,^ and the A 3 B 5 
compound bulk inversion asymmetry (BIA) mechanism described by the Dresselhaus term.— 
In the principle crystal axes, the bulk Dresselhaus 3D spin-orbit interaction Ti^, can be 
written as^ 

Hd = S(cr ■ k), where k x = k y k x k y — k z k x k z , (2) 

k y and k z can be derived from k x by cyclic permutations, and 5 is a parameter. Here 
kj (j — x,y,z) are the projections of the momentum operator k = — i V + eA /he of an 
electron, A is the vector-potential of the magnetic field B (9, <p) , 9 and <p are the polar 
angle and the azimuth of B , and — e is the electron charge. In a strong confinement limit, 
when carriers are in a quasi-2D regime, Ti.£> reduces to a 2D Dresselhaus Hamiltonian. For 
a rectangular [001] QW of the width d , 

H D = oi D (a x k x - (Tyk y ) (3) 
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with an = —b(k 2 z ) = —5{n/d) 2 . 11 i 12 i 13 i 14 i 15 The 2D Dresselhaus Hamiltonian possesses a 
discrete symmetry C 2v of the [001] face of a cubic lattice of the T d symmetry typical of 
the zinc blende modification of A 3 B 5 crystals, while the Rashba SO Hamiltonian 

H R = a R (a x k y - a y k x ) (4) 

possesses a continuous group C ^ of the rotations about the [001] axis. Experimental data 

nnq 

suggest that for GaAs QWs both SIA and BIA terms are of the same order 
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of magnitude. It is a general consensus that in narrow-gap compounds the SIA mechanism 
typically dominates,— however, the ratio of the coupling constants of only a^/ao ~ 2 has 
been recently reported for InAs QWs of width d — 15 nm.— The tunability of by gate 
voltae believed to be of critical importance for the operation of semiconductor 

spintronic devices. 

The spatial dependence of Zeeman energy Hz(r , <x) stems either from the inhomo- 
geneity of the field B — or from a position-dependent g -factor . 2 ^ 2 ^ 2 ^ 2 ? Recently Kato et 
a/.— achieved operating electron spins in parabolic Al x Ga i_ x As QWs through the g -tensor 
modulation technique based on the difference in the spatial dependences of the various com- 
ponents of g = g(z) . Graded parabolic QWs were originally designed for producing a high 
mobility electron ga o 30 i 31 by applying the modulation doping technique.— The early work on 
their dynamical responses was mostly concerned with the effect of electron-electron interac- 
tion on the electron orbital dynamics and the related transition frequencies,^ and also with 
generalizing the Kohn theorem.— The recent success in electrical spin manipulation^ shifted 
the interest to the spin flip transitions in such systems and the effect of various mechanisms 
of SO interaction. 

Thd existence of several mechanisms of SO coupling makes it important to develop reli- 
able experimental techniques for identifying them. Also, the relative efficiency of different 
SO coupling mechanisms strongly depends on the choice of the semiconductor materials and 
the shape of a QW or a heterojunction. The efficiency of (7 -tensor modulation technique 2 ^ 
is based on the anomalously small g -factor value in Al x Gai_ x As, \g\ < 0.1 , and is there- 
fore specific for GaAs based devices. Developing similar techniques for narrow-gap A 3 B 5 
semiconductors with typically large g values, \g\ ^ 10, needs different approaches, and we 
show that orbital mechanisms of SO coupling can be rather efficient for them. Moreover, 
we prove that the efficiency of different mechanisms depends strongly on the polarization of 
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the ac electric field E . 

To consider all these problems in a framework of an unified approach and to derive 
analytical expressions for the transition probabilities, we make two basic assumptions. First, 
we choose a parabolic shape for the QW. Second, we accept that SO coupling is not too 
strong, i.e., the SO coupling energy E$o is small compared with all different energies, 
including the confinement energy hu , the cyclotron energy H(J C , and the Zeeman energy 
hu s . Here ujq is the characteristic frequency of the parabolic potential, u c = eB/mc is the 
cyclotron frequency of electrons with the effective mass m, and uj s = gfiBB/h is the spin 
resonance frequency, with [is = ehjlm^c being the Bohr magneton. These assumptions 
allow one to account for Hso only in the matrix elements of spin transitions and disregard 
the effect of Hso on the position of energy levels.— In this approximation, one can take 
advantage of the exact solutio n 38 ! 39 for the quantization of electron levels in a magnetic field 
tilted with respect to the confinement plane. 

Meanwhile, applying the exact solution found by Maan^ 8 - and Merlin^ 9 , is not straight- 
forward when it comes to calculating matrix elements of spin flip transitions. Indeed, their 
solution depends on an auxiliary angle 7 defined by the decoupling condition of two normal 
modes, hence, all matrix elements depend on 7 . We have found an extensive set of remark- 
able identities relating 7 , frequencies oj^iO) of two eigenmodes, and the polar angle 9 . 
The symmetry of the problem underlying these identities is far from obvious, but the identi- 
ties permitted us to eliminate 7 and find all final results in an explicit form. Some of these 
results, without the derivation, were published in our previous papers . 40 ' 41 We expect that 
the technique will facilitate developing the theory of different properties of parabolic QWs 
which is timely because of the recent progress in experimental wor k 25 ' 29 and the competition 
of the different mechanisms of SO coupling that manifests itself in various phenomena. 

In this paper we develop a general theory of EDSR in QWs which is caused by the 
standard Hamiltonians of the H or i,( k , a ) type for two basic geometries: with the ac electric 
field in the QW plane and perpendicular to this plane. We solve exactly the problem of an 
electron confined in a parabolic QW being a subject to a tilted magnetic field B and find 
the EDSR intensity in the Dresselhaus and Rashba models vs the B direction. Our results 
show that electric dipole spin resonance is especially strong when it is excited by an in-plane 
electrical field. However, we show that it is also strong enough in the geometry when the 
time-dependent potential is applied to a gate. Our results demonstrate convincingly that an 
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efficient electrical spin manipulation can be achieved through the orbital mechanisms of spin- 
orbit coupling. Our results also suggest that the angular dependence of EDSR intensity is an 
unique characteristic of various competing mechanisms of spin-orbit coupling contributing 
to EDSR. 

The paper is organized as follows. In Section II we develop a theoretical approach to an 
electron in a parabolic quantum well and derive operator expressions for the basic variables 
like the coordinates and kinetic momenta. While we apply these results for calculating the 
EDSR intensity, they are rather general and can be applied to different problems related to 
parabolic quantum wells. In Section III we calculate the EDSR intensity for the Dresselhaus 
and Rashba 2D spin orbit coupling Hamiltonians for the two basic geometries with an in- 
plane and perpendicular-to-plane electric field. We also investigate in detail the dependence 
of the EDSR intensity on the direction of the magnetic field and present figures that illustrate 
the basic mechanisms and characteristic features that should help assigning specific EDSR 
bands when experimental work begins. In Section IV we develop a theory of EDSR in a 
parabolic well with a 3D Dresselhaus SO interaction and unveil the specific features of EDSR 
that distinguish it from the results of the 2D model corresponding to a strong confinement 
limit. The discussion of obtained results and estimates of the EDSR intensity are given in 
Sect. V. Appendix A includes a number of identities that have been highly instrumental 
in deriving the results of Sections III and IV, and that we expect to be useful for future 
theoretical work on parabolic quantum wells. Appendix B supports the calculations of Sec. 
IV. 

II. AN ELECTRON IN A PARABOLIC WELL A SUBJECT TO A TILTED MAG- 
NETIC FIELD 

Everywhere in what follows we suppose that SO corrections to the energy levels are small 
as compared with the separation between adjacent Landau levels and Zeeman sub levels. 
Therefore, we disregard the effect of SO coupling on the energy spectrum and begin with 
the Hamiltonian of the orbital motion in a parabolic quantum well 
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where A is the vector-potential of the magnetic field B . Because the Hamiltonian is 
quadratic in the coordinates, the problem can be solved exactly, and the energy spectrum is 
well known. 38 39 However, because for calculating matrix elements of EDSR we need explicit 
expressions for the operators of coordinates and momenta, somewhat lengthly calculations 
should be performed. 



A. Diagonalization of the Hamiltonian 



First, we perform a transformation U from the original coordinate system x, y, z related 
to the principal crystal axes to a new ("primed") reference system x', y', z' with the z' axis 
parallel to B and the y' axis lying in the x, y plane. The coordinates in both systems are 
related as: 

f x'\ 

where 



y 



u 



u 



y 

\ z ' J 

cos 9 cos — sin sin 9 cos <t> ^ 
cos 9 sin <fr cos (f) sin 9 sin 
— sin 9 cos 9 



(6) 



Similar equations are valid for momenta projections, fcj = Uu'ki' , where i = x,y,z and 
summation over repeated coordinate indices is implied. Because y' does not appear in the 
potential energy mujQZ 2 /2 , it is convenient to choose the Landau gauge with A x i = 0, A y i = 
Bx', A z i = . Then y' is a cyclic variable, the Landau momentum k = k y > is a c -number, 
and 



Hn 



= ^[-di-di + ik+x'/xr] 

+ -mco>Q(— sin^ x + cos 9 z') 2 , 



(7) 



where A = ych/eB is the magnetic length. To decouple the motion in two degrees of 
freedom in Eq. (jJJ), we perform a rotation 



x' = cos 7 £ + sin 7 rj, 



sin 7 £ + cos 7 rj. 



(8) 
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Then the cancelation condition for the mixed product £rj results in an equation on the 
auxiliary angle 7 

sin 2 7 = {u /uj c ) 2 sm[2{e + 1 )}, (9) 

and the coefficients at £ 2 and rf , together with the kinetic energy, define the eigenfrequen- 
cies 

u;|(6>) = cj 2 cos 2 7 + u; 2 1 sin 2 (6 l + 7) 

uo 2 (6) = ^ 2 sin 2 7 + cc> 2 cos 2 (# + 7) (10) 

of two normal modes, £ and rj . Equations and (jlOj) complete the diagonalization of the 
quadratic part of Hq and are in agreement with the results by Maan^ and Merlin.— 

Eliminating in Eq. (J7J) the term that is linear in x' can be achieved by shifting £ and 77 

by 

— kX 2 cos 7 



£0 



cos 2 7 + (loq I uj c ) 2 sin (6 + 7) 

= -\ 2 k C ° S(g + 7) , 
cose* 

—k\ 2 sin 7 
^° sin 2 7 + (c<j /u; c ) 2 cos 2 (6* + 7) 

= -X 2 k S J^±J1. (ii) 
cos# 

In these equations, first expressions for £ and 77 come directly from calculations, while 
the simplified form of them can be found by using Eqs. (jA13|) . 
The free term in Eq. (J7j) equals 

Mk) = |^ - y + M) = 0. (12) 

Vanishing of this term follows from Eqs. Q - (fTTj) . and independence of the energy from the 
Landau momentum k is a consequence of the translational symmetry in the (x, y) plane. 
The shifts £0,770 satisfy the identities 

cos 7 £ + sin 7 ?7 = -X 2 k, 
sm{6 + 7) £0 - cos(# + 7) Vo = . (13) 

Equation (fTUj) allows one to find the spectrum of the system, i.e., the set of two frequencies, 
{u^(9), uini.0)} ■ However, because Eq. (fTU|) includes the auxiliary angle 7 that should be 



found from Eq. the shape of the curves oj\{9) and ui 2 {9) is not obvious, depends on 
the ratio uj c /u)q , and identification of a specific eigenvalue as oj 2 {9) or uj 2 {9) is a matter of 
convention. In what follows, we identify ct>$ (0) and uj v {9) as the frequencies that at 9 = 
coincide with the cyclotron frequency uj c and the confinement frequency u> v , respectively. 
Then, using Eqs. (|A14|) and (jA16|) . we derive explicit expressions for oj 2 {9) and uj 2 (9) : 

w£(0) = ^ (u)l + lo 2 c - fi 2 sign{u; - ^ c }) , 

u 2 {9) = X - (u 2 + uj 2 + fi 2 sign{^ - w c }) , (14) 

where 

ft 2 (#) = vW + ^e -2cj 2 w c 2 cos20. (15) 

This spectrum is plotted in Fig. QJi as a function of cj c / c^o for three values of the polar 
angle 9 . For 9 = 0, the spectrum consists of the cyclotron and confinement branches, 
£ and rj , respectively. However, for 9^0 the branches interchange at u c = ujq ■ Left 
parts of the lower branches and right parts of the upper branches are described by the 
solution and show a cyclotron-like behavior. The opposite parts of the same branches are 
described by the u v solution and show a confinement-like behavior. Of course, all branches 
are continuous and smooth. 

The same spectrum is plotted in Fig.^D as the function of 9 for three values of oj c /ujq . All 
branches retain their identity, £ or rj , in the whole interval < 9 < ir/2. All rj branches 
originate at the same frequency ui n {9 = 0) = u ; they are upper branches for u c < ojq and 
lower branches for u c > uj . On the contrary, £ branches are lower branches for u c < uj 
and upper branches for uj c > u ; they originate at the frequencies uj c that are B -dependent. 
For <jj c = to , the spectrum is described by a simple equation to±(9) = cu (l±sm 9) 1 / 2 . These 
curves are separatrices dividing £ and rj regions, hence, neither £ nor rj identity can be 
ascribed to them. 

The Hamiltonian H of Eq. ((7j) , when written in the variables £ and rj , reads 

^o = E[-^? + ^c(C-Co) 2 ], (16) 

where the summation is performed over £ = £, r\ . After introducing step operators 

C - Co = \fh/2mu ( (a^ + a c ) , 

kc_ = —id^ = i\J~mwj2h{a^ — a^) , (17) 
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the Hamiltonian takes the standard oscillator form 

= 7;^2^ UJ d a t a ( + a c a t)- ( 18 ) 
1 c 

B. Operator representation for coordinates and momenta 

Similarly to Eq. (pj). one can express the components kj = — zV 'j + (e/Kc)Aj of the kinetic 
momentum in the crystal frame through its components in the primed frame. In the latter 
frame, the vector-potential A contributes only to the component &y that equals 

k y , = k + x'/\ 2 = [cos 7 (f - Co) + sin 7 {rj - r] )]/\ 2 . (19) 

When deriving Eq. (JTHJ), the upper identity of Eq. (JT3j) has been used. Using Eq. (jHJ), one can 

express k x i and k z < through the derivatives and rewrite the components of the momenta 
in the crystal frame in terms of a^, by applying Eq. (|17|) 

k x = J2( X t a < + X C a t)' 
C 

c 

k z = ]T(Z c a c + Z*a+). (20) 

C 

Here 



*e = 


V 2ft 




-z cos( 


9 + 7) COS (f 


) cos 7 sm 


X v = 


lmuj v 
V 2ft 




-z sin(f 


? + 7) cos 


) sm 7 sm 


n = 


V 2ft 


z cos( 


9 + 7) sin </ 


> cos 7 cos 




jmu! v 
V 2ft 


z sin(f 


? + 7) sin$ 


U c ■ ± 

sm 7 cos (p 



Z, = i^sin( S + 7 ), 

z - = _ Vilr cos(e+T) - < 21 > 

These coefficients depend only on the angles (8, <ft) and the eigenfrequencies uj^ and do 
not depend on k , £ , and r] . Therefore, despite the fact that Eqs. (J20|) and pTf have 
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been derived in the Landau representation, the final results are gauge invariant and may be 
conveniently employed for finding operators (x, y, z) of the coordinates (x, y, z) . 

It is seen from Eq. (J5J) that the operator z can be found in terms of (a^, a£) from the 
difference H — h 2 k /2m. Employing Eqs. (|18|) . (|2(J|) . and (|21j) and some of the identities 
of Appendix A, we arrive after somewhat lengthly algebra at 

z = ^(ih/rruj < ;)Z c (a c + a{). (22) 
C 

The sign of z has been found from the commutation relation [k z ,z] = —i, and one can 
easily check that z = (i/h)[Ho, z] = hk z /m. 

For finding x and y , we generalize the Johnson and Lippman 42 procedure for the quan- 
tization of electron motion in a strong magnetic field for an electron confined in a parabolic 
quantum well. Using the equations of motion for the operators of momenta and coordinates 

k x = ^c(byk z b z ky), ky = u) c (b z k x b x k z ), (23) 

x = hk x /m, y = hky/m, z = hk z /m, (24) 

where b = B /B , and eliminating the components of the momentum k , we arrive at the 
equations 

± = (X 2 l y + bj)/b z , y = (-X 2 L + b y z)/b z . (25) 
Equations (|2l)|) suggest existence of two integrals of motion 

x = x - {b x z + \ 2 k y )/b z , y = y- {b y z - \ 2 k x )/b z , (26) 

that generalize the well known guiding center coordinates. They obey the same commutation 
relation 

[x ,y ] = i\ 2 /cos6 (27) 

and commute with operators . Remarkably, in a parabolic well the operators xq and 
yo involve the vertical coordinate z . Nevertheless, Eq. (127)1 ensures that the degeneracy 
of states is completely controlled by B z , the component of B perpendicular to the QW 
plane, and is described by the standard Landau formula n^O) = cos#/27rA 2 . 

Because the operators k x , k y , and z are already known, Eq. (j2Hj) allows one to find x 
and y . To derive explicit expressions for these operators, it is convenient to use identities 

+ {iuj c / u)() sin 9 cos <pZ^ = (iu c / u^) cos 9X^ , 

+ (iuic/u^) sin#sin0Zf = (iu c /u^) cos 9Y^ (28) 
11 



that can be checked using Eqs. (fTUjl and (}2*Tj) . Finally, the operators of in-plane coordinates 
are 



c 



(29) 



C 



Because Zj- = —Z^ , Eqs. f!29|) are similar to Eq. (|2*2*j) and differ from it only by presence of 
the guiding center operators. 

Therefore, six operators of coordinates and momenta, Xj and kj , j = x,y, z , are ex- 
pressed in terms of four Bose operators (a^,a^) and two generalized guiding center coor- 
dinates (x , y ) . The latter commute with all operators (a^, a^) ; hence, their presence in 
Eq. (J29|) does not influence electron dynamics in homogeneous external fields that are cou- 
pled to an electron only through the operators (a^, a~£) . Operators x j and kj , defined by 
Eqs. (J2IIJ), (J22), and (J2HJ), obey the standard commutation relations 



where ej£ m is the Levi-Civita tensor. 

III. INTENSITY OF EDSR 

The EDSR occurs in QWs due spin-orbit Hamiltonians Hso described by Eqs. (J2J) - 
(0J) which mix the electron spin projections on the magnetic field direction, hence, each 
electron state acquires admixture of the opposite spin projection that is small when Hso 
is weak. As a result, a time dependent electric field E(t) causes EDSR, i.e., electrically- 
induced spin-flip transitions. Matrix elements of EDSR are comprised of two contributions, 
of which one comes from the perturbation of the wave functions of stationary states and 
the second from the direct coupling of the electron spin to the field E(t) . A convenient 
way for calculating EDSR matrix elements is based on employing a canonical transformation 
exp(T) eliminating H so in the first order of perturbation theory.^ After the transformation, 
the time-independent part of H conserves the electron spin projection on the magnetic 
field direction, and spin-flip transitions are induced only by the spin-orbit contribution 
e(v so ■ A)/c to the time-dependent part of H , where A (t) is a time-dependent vector 
potential, and v so is the spin-orbit part of the velocity operator. 



[kj,Xi] = -i5j£, [kj,k e } = -i\ 2 e j£m B m /B 



(30) 
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Because of the equation of motion v = (i/K)[H, r] , the matrix elements of the oper- 
ators v and r for spin-flip transitions are related as (v) = iu s ( r } , with u s being the 
spin-flip transition frequency. Neither of the two competing contributions to v includes 
the factor u s , hence, it indicates the existence of massive cancellations that tremendously 
complicate calculations based on the operator v . It is much more convenient to write the 
time dependent part of the Hamiltonian as H int (t) — e( r ■ E (i)) . In the original represen- 
tation, the coordinate operator r is diagonal in spin indices, and H int (t) produces spin-flip 
transitions due to the level mixing. However, after the T -transformation the operator r 
acquires a SO part f so = [T, r] that drives spin-flip transitions. 

The total Hamiltonian of an electron confined in a parabolic QW is H = H + Hz + 
H so + H int {t) , with Hz = \q^b{, cr ■ B) . The energy levels of the Hamiltonian H + Hz 
are 

E a (nz,nr,) = Y,toc(6)(n c + l/2) + hu; 8 (T/2, (31) 
C 

where tV > and the spin index a = ±1 . The spin-flip frequency uj s should be taken 
algebraically; uo s < for electrons in negative g -factor semiconductors. The time indepen- 
dent spin-orbit interaction H so will be considered as a perturbation. 

The term H so — Hd, Hp, or Hr in the time-independent part of the Hamiltonian H 
leads to the mixing of spin sublevels. As we have stated above, it can be eliminated in the 
first order of the perturbation theory by a canonical transformation exp(T) . The operator 
T is nondiagonal in the orbital quantum numbers (n^,n v ) , and its matrix elements are 

(n g ,n a^n^n^a) = — s / — -. (32) 

In terms of T, the matrix elements of spin-flip transitions diagonal in (n^,n v ) are 

(n^n^T \(E ■r BO )\n^n v ,i) 
= J2 {(n^n V} ] \ f\n' v n' n ,i)(n'^n' n \{E-r)\n^n v ) 

- (n^n v \(E-r)\n'^n' v )(n^,n' v ,^ \f\n^n v ,i)} , (33) 

with n'^ = ± 1 and n' v = n v ± 1 . The sum in Eq. (|3*3*|) is restricted because the operator 
r is linear in (a^, a£) , cf. Eqs. (|2*2*|) and (|29jl. This restriction significantly simplifies matrix 
elements (n'^n' , j {Tln^n^ j) since their denominators do not depend on (n^,n v ) and are 
equal to h(±uj^ + u s ) . 
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When the operator H^o is linear in momenta and consequently in the operators (a^, a£) , 
Eq. (J3~3*)) reduces to 



+ 



(n f ,n,,t \{E-r so ) \n^,n v ,i) 
_ _ UE-l* ( )(0 o UH S o\l c ,i) 

(E-l ( )(l o UH so \0 o i) 



with 



ih ih 



ih 



(34) 



(35) 



Remarkably, this matrix element does not depend on (rig, n v ) . For calculating the matrix 
element (| \H so \ |) in Eq. ()34jl . we transform Pauli matrices from the crystal frame (x, y, z) 
to the primed frame (x',y',z') as &i = , similarly to Eq. (|B|). Substituting cr^ into 

the 2D Dresselhaus and Rashba Hamiltonians of Eqs. (J3J) and (|4*jl. and using (j lovl |) = 1 , 
(T l<vl I) = ~~ * i an d (T I) = , we arrive at 



(T |4o| I) = E(^c + ^c +a c) > 



(36) 



with the matrix elements fZ^ that should be found from Eqs. (|20|) and (|21|). Applying 
Eq. ()A13|) . we find for a Rashba Hamiltonian 



H R,C = OLR^jmuj^llh {u c /u c ±l)f ( , 
with /e = cos(8 + 7) and f v = sin(6 l + 7) , and for a 2D Dresselhaus Hamiltonian 



(37) 



H. 



2D,C _ a D 



1 cos. 



— sin 



2H 



uv cos # 



=F cos# 



fe- 



rn 



Substituting Eq. (jHEJ) into Eq. (pHjl . we find 



x H[# c -(z c -£r + tf+(z c -£)] 

+ ^i?^- •£)]} 



(39) 
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A. Electric field perpendicular to the QW plane 



Let us start with a time-dependent electric field perpendicular to the QW plane, 
E (t) || z . In this geometry = 1*^ , and the matrix element of Eq. (|39|h L z = (n^, n v , | 
Izsoln^rirj, I) , equals 



+ w s 



s 2 9 



For a Rashba SO coupling, substituting Eq. (|37j) into Eq. (j4*Uj) results in 



L^4l + - 

R 2n v w, 



sin 2(0 + 7) 



2>(0) 



with 



(40) 



(41) 



V{6) 



[OJc 



2 -uj 2 ) 

r\ s J 

It, .2 



2 2 2 n 2/2,2 2\ 

uj u c cos 9 - uj s {uj +uj c - u s ); 



(42) 



here, Eq. ()A10jl has been applied. By using Eq. (|A15jl . one arrives at the final resulb^ 1 ^ 

a R uj s (lo c + u s )u c 



sin 29 



(43) 



2h V{6) 

For Dresselhaus SO coupling, a similar procedure leads to the matrix element of Ref. j^J 



J D 



sm 20 cos (a; c — o; s J 



ft £>(0) 

— z cos 20 (u c cos 2 6 — uo. 



(44) 



B. Electric field in the quantum well plane 



For an in-plane electric field, we calculte the scalar product (I ^ ■ E) . If the field 
E (t) is polarized at the angle ip to the crystal x-axis, E (t) = E (t) (cos i/j, sin ip, 0) , then 
(l c ■ E) = lfE(t) with 



h 



\ 2mo; ( 



cos(0 — ip) — i — - sin(0 — ib) 

lu c cos 



(45) 
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When deriving Eq. (JUJ, the identity of Eq. (jA7|) has been used. Equation (}4T)j) allows one 
to rewrite the matrix element = (n^,n v ,^\ \x so cos ip + y so sin -0| rig, n v , J.) of Eq. as 



+ EM?) 



uj}-uj 2 s 



u 2 - u 2 



g£ - g 

W C 2 2~~ 



(46) 



For Rashba SO coupling, using Eq. (pTTj) results in 



4 



— I- 



a.R cos(0 — ip) 
% 

. a R sin(0 - ip) 

COS ^/idX, 



E 

c 

1 + 



1 /<v 



u c u 2 - u 2 



+ to. 



u 2 - u 2 



10, 



f 2 , 1 2 

c u}-u 2 



E 



(47) 



Applying identities of Eq. (|A3J) . one finds an explicit expression for the dependence of the 
matrix element on the magnetic field direction 

L R = -^COS(0-V) 

uj c cos 2 8(uJq — uj 2 ) + uj s (lJq + u 2 sin 2 6 — uj 2 ) 

x W) 

(u c + u s )(ujI - a; 2 ) 
V{9) 



— i-^- cos6 sin(0 — ip) 



(48) 



In the strong 2D limit, when uj ^> uj c , this expression simplifies significantly, and we arrive 
at the result of Ref. |4l] 



JR ~ 

+ i sin(0 — ip) (uj* + cj s cos 



cos(0 — ip) (ui* cos 6 + u s ) 



(49) 



Here uj* = uj c cos # is the cyclotron frequency in a tilted magnetic field; from now on 9 is 
defined as < 6 < vr/2. 

For a 2D Desselhaus SO coupling, substituting Eq. (|HK|) into Eq. (f4T)J) results in 



L%/a D = 
cos(0 — -0) 

h 



E 

c 



cc^ i cos 20 — sin 20 cos # 



ijj r COS 6 1 



+ 



u% - UJ 2 



,sm. 



I COS COS I 
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isin(0 — VO ^C^C u; s («cos20 — sin 20 cost 

2^ ,.,2 _ ,,,2 



cos 9hu c ^ oof — ioi 



ujc cos 9 



+ sin 20 — i cos 20 cos #] . (50) 
By using identities of Eq. ()A3jl . an explicit expression for this matrix element can be found 



^ COs(0-^) 

L D /a D = — 



uj c cos9(lUq — uo s ) 

W) 



x (i cos 20 — sm 20 cos 9) 

{ujI + u; 2 sin 2 - w 2 ) . . 
+ cj s — ^— — — (sm20 - zcos20cos0) 



V(9) 

zsin(0-<0) (wg -^ 2 ) 



Wo Ucos z0 — sm 10 cos( 



ft V{6) 

+ u; c cos 6* (sin 20 — z cos 20 cos 6 1 )] . (51) 

In the strong 2D limit, ujq ^> uj c , through a number of transformations this expression 
can be brought to the final form^* 

L d = % , ,2 J — 2r[sin(0 + ^)K*cos0-^ s ) 

— i cos(0 + ip){oj* — uj s cos 0)] . (52) 

In conclusion, we note that in this and the previous subsection we considered the Hamil- 
tonians Hd and Hr on a similar footing. However, there is a considerable difference in 
their ranges of applicability. The Hamiltonian Hd originates from the 3D Hamiltonian H.D 
and is valid only in the strong confinement limit when d <C kp , A , where hp is the Fermi 
momentum. When the confinement length d becomes comparable to k F l or the magnetic 
length A , additional terms like k x k 2 y should be included in the 2D Hamiltonian. From 
this standpoint, only Eq. (}52|) has a real physical meaning while Eq. (151)1 should be only 
considered as an auxiliary one. Having those limitations in mind, Eq. ()51)1 will be applied 
in Sec. IIII C 21 only to electrons confined at the ground level, n^,n v = 0. Eq. (jH)) is a 
subject to similar restrictions, and EDSR with n^,n v ^ electrons in a field E || z will 
be considered in Sec. IIVI in the framework of a more general theory. 

The applicability range of the Hamiltonian Hr is much wider because it is applicable 
not only to zinc blende crystals under the conditions of 2D confinement but is inherent 
in the wurtzite modification of A 3 B 5 compounds as a bulk property.^ Currently, active 
experimental work on SO properties of microstructures including the wurtzite modifications 
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of InAs,— InN,— GaN , 46 ' 47 i 48 etc. is under way, and large SO splittings up to 9 meV have 
been reported.— Following the initial calculations of the band structure of wurtzite type 
compounds,— more general models have been developed recently , 50 ' 51 and it looks like the 
band folding in the hexagonal direction (resulting from the different size of the elementary 
cells in the zinc blende and wurtzite lattices) plays a role in developing large spin splittings.— 
Therefore, the Hamiltonian Hr will be applied for arbitrary n^,n v ^ 0. 

C. Polarization dependence of EDSR intensity: General properties 

In this section, we discuss at a qualitative level the basic properties of EDSR following 
from the equations of Sections IIII Al and IIII Bl 

In all cases, the intensity of EDSR shows a pole when one of the eigenfrequencies coincides 
with |u; s | , W({O po i e ) = \u s \ . It is seen from Fig. [T]d that such a pole always exists because 
one of the frequencies u^(0) vanishes for 9 = tt/2M For u > u c it happens to the ui^ 
mode while for ujq < uj c to the u v mode. The magnitude of ^poie depends on the relative 
magnitudes of u , u> c , and uj s . For u ^> u c , pole depends mostly on the ratio uj s /uj c . 
It is very close to 7r/2 for GaAs because of its anomalously small g -factor, and equals 
#poie ~ 0.89 x 7r/2 for InAs and 9 po \ e ~ 0.78 x tt/2 for InSb. Therefore, for the materials like 
InAs and InSb the maximum can still be achieved in the range of 9 and B values where the 
transformation of Eq. ({HZ} remains justified. Additional dependence of the matrix elements 
on 9 comes from the numerators of Eqs. (|43|) and f)44j) and their analogs for an in-plane 
electric field. This dependence is SO coupling specific. Singularities at the poles are cut-off 
by a level width. 

The 9 dependences of the EDSR transition matrix element for J57 || z and E _L z are 
significantly different. First, for the perpendicular-to-plane polarization, the factor sin26 l 
in Eq. (|4"H|) and the factor sin 9 in Eq. ()44|) nullify these matrix elements at 9 = 0. There 
are no such factors in the matrix elements of Sec. IIII Bl Therefore, EDSR can be observed 
only in a tilted magnetic field B when E (t) \\ z , but an arbitrary orientation of B can 
be used (including B || z ) for an in-plane E (t) . 

The second difference in EDSR with different E (t) polarizations concerns its intensity. 
This difference is essentially pronounced in the strong confinement regime, ujq ^> uj s ,uj c . 
Comparing Eqs. (|4*H|) and (|4*4*|) with Eqs. (|4"H|) and (|52"|) . one can estimate the corresponding 
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matrix elements as ~ a^o^s/fi^o ^ OT E(t) || z and as ~ a^o/^c for E (t) _L z, 
respectively, when the inequality u c > \u s \ is satisfied. The first estimate indicates that 
when E (t) || z , EDSR is possible only due to the deviation of the system from the strict 
2D limit. On the contrary, when E (t) 1 z, EDSR survives in the strict 2D limit. The 
ratio of the matrix elements is about uj c uj s /uJq -C 1 , hence, in the strong confinement limit 
an in-plane field E (t) is much more efficient than a perpendicular-to-plane field. 

For SIA, it is also instructive to compare the magnitudes of the EDSR transition matrix 
elements in a strong 2D confinement limit, L% , described by Eq. ()52j) . and in the bulk, 
L3D ~ 5m/h 2 , as found by Rashba and Sheka£. We find L^/L ~ ujq/u c ^> 1 using 
an ~ 5mu) /h from Eq. ()67|) below. This enhancement of EDSR in a QW in in-plane 
geometry can be attributed to a strong confinement with the square of the confinement 
momentum, muo/h, large compared with A~ 2 . Extraordinarily high efficiency of EDSR in 
the in-plane geometry has been emphasized in Ref. 0]. 

We note that despite the fact that the intensity of EDSR with E (t) || z is much less 
than with E (t) _L z , it is usually high enough for efficient electrically manipulating electron 
spins. Specific estimates can be found in Ref. and will be given in Sec. IVl below. 

As has been discussed above, there is a similarity in some of the properties of EDSR 
caused by BIA and SIA. The most striking difference in the effect of these mechanisms is 
seen in the angular dependences of the EDSR intensity, especially in its dependence on the 
azimuth (ft , as will be discussed in Sec. IIII C 21 



1. Dependence of EDSR on the polar angle 

Using the equations of Sec. IIII A| we will provide and discuss here the dependence of the 
EDSR intensity on the polar angle 9 as applied to InAs. We restrict ourselves to SIA because 
Eq. (jSJ) for 2D Dresselhaus coupling is applicable only under strong confinement conditions 
when ojq 3> uj c ; cf. Sec. IIII Bl A theory of BIA controlled EDSR will be discussed in Sec. IIVI 
in the framework of a more general approach. 

It is seen from Fig. [I]p that the energy spectrum is rather different for u c < ujq and 
uj c > ujq . Therefore, these two cases will be considered separately. 

The angular dependence of the EDSR intensity is controlled by the dependence of the 
square of the matrix element L Z R on 9 , Eq. ()43|) . and by the population difference of the two 
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spin sublevels that depends on the filling factor v{9) = n e /ni(9) , n e being the concentration 
of 2D electrons. To cut-off L Z R (6) near its pole and to find a realistic estimate of EDSR 
intensity vs 9 , we introduce a phenomenological level width T . 

When uj c ^> u s , one can consider a single spin-split level E ± hu s /2 with < v < 2 . 
For a Lorentzian level shape, the difference Au of the filling factors of two sublevels equals 



7T 



aIctm (ir ;3 «) + W2 



— arctan 



r 

(77 - Eq) - fao s /2 



r 

where r](i>) is the chemical potential that can be found from the equation 

V(v)-E _ 1 



(53) 



tan 7r(V — r 



+ signjz/ — 1} 



x 



1 fhuj 



(54) 



\ sin 2 7r(z/-l) V 2r 

The dependence of Av on v is shown in the Fig. El For small T , this is a triangle. 
Its vertices can be considered as cusps of the Au vs v curve. An up-cusp appears at 
v = 1 , and two down-cusps at v — and 2. Note, the up-cusp appears at odd v , while 
down-cusps at even v values. With increasing T , the cusps are smeared, however, they 
manifest themselves as a well pronounced maximum and two minima even for T as large as 
T = Hu s /2 . Below, we will find similar patters in the spectra of multi-level systems that are 
well pronounced even for a strongly tilted field B when u c cos 9 ~ u s . However, because 
of level crossings the regular alternation of up- and down cusps (and their correspondence 
to odd and even v factors, respectively) can be violated; an example of such an anomaly 
can be seen in Fig. |2b- 

In what follows, we use a Gaussian level shape exp[— (E — £'o) 2 /r 2 ] because it provides 
pronounced features in the intensity distribution that can be reliably assigned. This proce- 
dure will be applied to the denominator T>{9) of Eq. (f4*3~j) where it cuts the pole and also 
when calculating the populations of different orbital levels and their spin sublevels. How- 
ever, we do not renormalize the frequencies u c and u s in the numerator of Eq. (J43|) . and a 
similar procedure will be applied everywhere below. 

In Fig. Ok is shown the dependence of EDSR intensity on 9 for oj c = uj /2 and the 
electron concentration n e = 2n^(6 l = 0) . Therefore, for 9 = both spin sublevels of the 

= 0,12^ = level are filled, and EDSR develops only when B becomes tilted, ul{9) 
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decreases, and the > l,n v = level is getting occupied. The general shape of the curve 
is dominated by a pronounced maximum that is achieved near the pole of the denominator, 
uj^ = \u s \ , and shifts to lower angles with increasing T . Strong suppression of the EDSR 
intensity for small angles and for 9 « it/ 2 originates from the factor sin 29 in Eq. (J43j) . 

Assigning specific features in Fig.Ek can be done by following the populations of successive 
levels and the behavior of the chemical potential 77 with increasing 9 ; 7](9) is shown in 
the insert. For small V values, populating level begins only after both spin components 
of the previous level, with the quantum number (n^ — 1) , are completely occupied. This 
regularity holds up to v = 10 and manifests itself in regular tooth-like pattern of the 
curve r)(v) for small level width, T/huo = 0.01. The same regularity can also be seen 
from the ^-dependence of the EDSR intensity for T/hujQ = 0.05. Similarly to Fig. |2l 
up-cusps correspond to odd filling factors v , while down-cusps to even v . The larger is 
the gap between the successive energy levels, the larger is the "jump" in , and the 
more pronounced is the corresponding cusp. Small shifts of the cusps from integer v- 
values originate from the final level width T . Because of the general growth of the EDSR 
intensity with 9 at the left slope of the principal maximum, up-cusps at this slope are shifted 
to the right from odd v -values while down-cusps are shifted to left from even //-values. 
Remarkably, all these features are distinctly seen in the intensity pattern for T/hujQ = 0.05 
despite the fact that the function r/(z/) is already rather smooth for this level width. 

We have calculated only transition intensities and did not calculate the corrections to the 
g -factor that originate from the band nonparabolicity and SO coupling. These corrections 
should depend on the quantum numbers n^,n v . Because of the regular level alternation 
described above, only a single spin-flip frequency should be seen inside each window confined 
between two successive even values of v ; values for specific windows are indicated in the 
insert to Fig. 3a. Two spin-flip frequencies can be seen simultaneously only inside narrow 
regions of 9 where v passes through even values. 

The EDSR intensity for uj c /ujq = 2 and the electron concentration n e = 2n^{9 = 0) is 
shown in Fig. Eb; in this case, u v is the lower spectrum branch. The basic shape of the 
spectrum is the same: it is dominated by a strong maximum at uj v ~ \u s \ . However, there 
is a considerable difference in the fine structure. First, because the principal maximum is 
achieved at a lesser value of 9 , fine structure is distinctly seen on both sides of the maximum. 
Second, level intersections arise at relatively small quantum numbers n v ; the intersection of 
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(1, 1) and (2, j) levels is the first one. However, of much more importance is the intersection 
of (2, J,) and (3, j) levels because Fermi energy rj passes through the intersection point. As 
a result, the simple regularity in populating successive levels does not hold any more, and 
alternation of up- and down-cusps is violated. Between the well pronounced v = 5 up-cusp 
and the v = 7 down-cusp there exists a smeared v = 6 up-cusp. For this reason, signs of 
the following cusps change: a strong odd- v cusp v = 7 turns into a down-cusp, while an 
even- v cusp v — 8 becomes an up-cusp. 

The intersection of (2, j) and (3, |) levels happens when 3^ — |c<j s |/2 = 2^ + \u s \/2, 
i.e., it coincides with the zero of the denominator of Eq. (|43jl and the principal peak of 
EDSR. Because this intersection point also coincides with the Fermi level, to the left from 
this point only the (2, f ) — > (2, j) transition to the half-populated (2, j.) level is allowed. On 
the contrary, to the right from the intersection point both the (2, f ) — > (2, j) transition to 
the empty (2, j) level and the (3, f ) — > (3, |) transition to the half-populated (3, 1) level are 
allowed, hence, the intensity of EDSR is expected to increase abruptly by a factor of three. 
The inflection point seen between the v — 5 cusp and the maximum of the curve reflects this 
discontinuity smeared by the level width T ; smearing is rather strong because the slopes of 
(3, |) and (2, 1) levels plotted vs 9 are very close. The asymmetry of the peak with respect 
to the LO n = \u s \ point, including shifting its maximum to the right from this point, is a 
different manifestation of the effect of the increase in the number of electronic channels at 
the intersection point. Also, the coexistence of two transition channels, (2, f ) — > (2, J.) and 
(3, |) — » (3, |) , manifests itself in the high EDSR intensity of the background over which 
the even- v cusp at v — 6 is hardly seen. 

The coexistence regions of different transitions are shown in the insert to Fig. Eb- Tran- 
sitions between spin sublevels of the levels n v = 2 and n v = 3 , coexisting in a wide region 
of filling factors, should manifest themselves through splitting the spin-flip line into two 
components with close g -factors. The same is valid for n. n = 3 and n v = 4 levels (for 
larger 9 values). 

All the above analysis is based on a phenomenological approach with a single level width 
parameter F . In the framework of such a theory, it describes both the widening of £ and 
7] energy levels and the width of the spin-flip transition lines. However, the real physics 
might be rather different. It has been shown by Mel'nikov and Rashba^ that impurity 
scattering results in a dramatic narrowing of spin-flip lines, and narrow EDSR lines have been 
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observed in dirty semiconductor materials by Bell long before. 56 The underlying mechanism 
is the dynamical narrowing that is also responsible for the D'yakonov-Perel' spin relaxation 
mechanism.— Therefore, when two or more wide £ or rj levels are partially populated, one 
should anticipate observing several narrow spin-flip lines in the EDSR spectrum. Indeed, 
two narrow spin-flip lines were observed in n-type inversion layers on InSb, one of them 
showing a Fano-type profile;— apparently, they were caused by EDSR, but the mechanism 
has not been specified. The detailed theory of such spectra should be model dependent. 

Therefore, EDSR spectra include rich and nontrivial information about energy spectra. 
Extracting it can be achieved by a detailed analysis of the spectra supported by calculating 
energy levels and their populations. 

2. Azimuth dependence of EDSR intensity 

In the previous section, we considered a perpendicular-to-plane electric field, E || z , 
concentrated on the details of the EDSR dependence on the polar angle 9 , and restricted 
ourselves with SIA. Because of the axial symmetry of the Hamiltonian Hr , the azimuth 
dependence of EDSR intensity was isotropic. 

In this section we provide the full angular dependence of EDSR intensity with a special 
emphasis on its azimuth dependence; and with this end in view we consider both basic polar- 
izations of the ac field E and both 2D Hamiltonians, Ho and Hr . Because the equations 
of Sec. IIII Bl derived for the Hamiltonian Hp are applicable only under the conditions of 
strong confinement, we accept that uj ^> uj c ,uj s , and that the electron concentration is low 
enough, hence, electrons populate only the lowest quantization level, n^,n v = 0, for all 
polar angles 9 of interest; this restriction will be removed in the more general theory of 
Sec. |TV| To suppress fine features of the 9 dependence, we choose a relatively large level 
width, T = 0.2c<Jo ; it suppresses the resonance at uo^ = \u s \ and enhances EDSR in the 
small 9 region. 

At first we keep E || z and consider the dependence originating from the Hamiltonian 
Ho . It possesses only the two-fold symmetry of the C 2 V group. However, the matrix 
element L Z D of Eq. (|^4j) is odd to tt/2 rotations. Therefore, the intensity of EDSR acquires 
the four-fold symmetry that manifests itself in Fig. However, the joint effect of Hp and 
Hr (the latter one is totally symmetric!) eliminates this additional symmetry. The effect 
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of the interference of BIA and SIA is seen in Fig. it is especially strong for «d = ±a_R . 
The degree of asymmetry should allow measuring the ratio old/cxr. 

Now we turn to the in-plane geometry, E _L z , where the symmetry of the indicatrix 
is additionally lowered because the azimuth ip of the field E establishes a new preferred 
direction in the (x,y) plane. Figures and are based on Eq. (|49|). Because of the 
isotropy of the Hamiltonian Hr , the EDSR intensity depends only on the difference of the 
azimuths ip and <fi of the electric field E and magnetic field B , respectively. Therefore, 
for (p = ip , the EDSR intensity does not depend on the azimuth as it is reflected in the 
rotational symmetry of Fig. EH- Remarkably, in this geometry Eq. (f4*H|) reduces to its first 
term that vanishes when io 2 c cos 2 9 + u s = . That is why the angular indicatrix of Fig. 
drawn for InAs with u s = — 0.17u c , consists of two sheets touching in a single point. The 
low- 8 feature became visible because of the large F = 0.2a;o . From this standpoint, it is 



instructive to compare the figures from our previous paper, Ref. |41|, with the set of figures of 
this section; they complement each other rather well. In Fig.Eb, the direction of electric field 
is fixed along the [110] crystallographic axis, and the <fi dependence of the EDSR intensity 
shows only the two-fold axis symmetry. 

Figures Et and 01 are based on Eq. (J52J) . For Fig. Et, E \\ B ± and therefore = ip , 
B ± being the projection of B onto the confinement plane. The matrix element of Eq. (|52|) 
has only a two-fold symmetry axis that corresponds to the symmetry of the Hamiltonian 
Hd ■ However, its square possesses four-fold symmetry that one easily recognizes in Fig. Efc. 
In Fig. Eli, direction of E is fixed as E || [110] , and the symmetry of EDSR intensity is 
reduced again to a two-fold axis. 

Figure illustrates the joint effect of Hr and Hd Hamiltonians for E \\ B±. This 
figure, drawn for = old , reflects the two-fold symmetry inherent in the BIA Hamiltonian 
H D . 



D. Experimental data 

EDSR is well documented in 3D where various mechanisms of it have been discovered and 
identified . 6l56l59l6 ° Experimental data regarding 2D systems are very scarce yet. Ironically, in 
2D the effect of spin-orbit coupling on spin resonance is best understood in Si/SiGe quantum 
wells where SO coupling is notoriously weak; 2 ?^ 1 the success is mostly based on long spin 
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coherence times. Electrical monitoring of electron spins in AlGaAs parabolic quantum wells 
near the g ~ point has been achieved by the (7 -tensor modulation technique. 29 It was only 
very recently that EDSR driven by the orbital mechanism has been reported for a A 3 B 5 
quantum well.— And ironically again, this observation has been made with an AlAs quantum 
well where the weakness of SO coupling manifested itself in a small g -shift, g ~ 1.99. 
Nevertheless, the intensity of EDSR exceeded the intensity of EPR (excited by the magnetic 
component of microwave field) by four orders of magnitude. As Schulte et al. emphasize in 
their paper,— the very fact of the observation of a spin flip line represented a puzzle because 
in their sample the EPR intensity was two orders of magnitude smaller than the noise level. 

We note that AlAs is not a propitious material for applying our theory because of the 
small effective mass, m/m = 0.46, and small g -factor, g = 1.99. As a result, even for 
magnetic fields B ~ 1 T, the cyclotron and spin frequencies are in the microwave range 
and the criterion of strong cyclotron quantization is at the verge of its applicability. In 
Ref. Idl, f c r p = 2.4 and f s r p = 1.1 for B = 1.219 T, and f c r p = 0.67 and f s r p = 0.31 for 
B = 0.3349 T, wher e f r s = u> c ^/2n and t p is the momentum relaxation time. Apparently 



for this reason, Ref. |§2j does not include any data regarding cyclotron absorption, and all 
the discussion is provided in terms of non-quantized orbital dynamics. The basic conclusions 
of Schulte et a/. 6 - are as follows. The signal of spin resonance originates from the effective 
magnetic field caused by the SO interaction. The polarization dependence of the resonance 
is well described by the effective field B c g = (2an/gnB){ k x z) following from the Hamil- 
tonian of Eq. (jH), and its intensity suggests ~ 5 x 10 -12 eV cm as a crude estimate for 
the SO coupling constant. 

The high intensity of EDSR in the B || z geometry, even when the sample was positioned 
close to the node of the microwave electric field E , was observed only when the field E 
was in-plane polarized; unfortunately, EDSR in a tilted magnetic field has not been studied. 
EDSR has not been seen with E || z These observations are in agreement with the 



predictions of Refs. and and the conclusions of Sec. |II 1 ( ' 

IV. 3D DRESSELHAUS HAMILTON! AN 

The physical parameter that allows reducing the Hamiltonian i-tp of Eq. (j2J) to its 2D 
form Ho is a small confinement length d or, what is the same, a large confinement fre- 
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quency ujq . Therefore, the expressions found in Sec. lIIII for the 2D Dresselhaus Hamiltonian 
of Eq. (J3J) and arbitrary values of lo c /ujq and (n^,n v ) are mostly of methodical interest. 
They have provided, side by side with equations found for the Rashba Hamiltonian of 
Eq. (J3J), an important outlook on the comparative strength of EDSR excited by in-plane 
and perpendicular-to-plane electric fields E (t) . However, a consistent description of 3D 
Dresselhaus systems can be achieved only by using the Hamiltonian Tio . Technically, cal- 
culating matrix elements of EDSR for that Hamiltonian is a challenging task. We will show 
in this section, as applied to a field E (t) || z , that using the operators of Sec. Ill Bl in 
conjunction with the identities of Appendix El is a powerful tool that allows solving the 
problem and deriving explicit expressions for matrix elements. 

In the primed reference frame, the term acquires a form Tip = 5<jj>Kf with ky = 
Ujj/kj , because the matrix U is orthogonal, hence, U^ 1 = U T . Of all the multitude of 
terms that enter in this expression, we need to select only those that contribute to the matrix 
element of Eq. (}3*3*j) . First, the selection rules in the spin operators show that 

(| \H D \ \) = 5K = 5{k' x , - ik' yl ) . (55) 

Second, because the operators of coordinates are linear in (a^, at) , the selection rules in 
the orbital operators in equations similar to Eq. select only those parts of the operators 
k'-, that change one of the quantum numbers, either or n v , by ±1 . Therefore, in what 
follows we denote by k'y these parts of the operators. Straightforward calculations in the 
crystal coordinate frame result in the following expressions for them 

C 

+ E (CW«C + C;.^a+) (2n c , + 1) , (56) 
CYC 

where = a^d^ . Coefficients C x -^' read as 

C x;CC , = Y C (X C ,Y$ + Y C X*,) - Z C (X C ,Z*, + Z C ,X*, ) 

+ X ( (\Y C \ 2 -\Z C \ 2 ) , (57) 

and coefficients C y .^i and C y .^i can be found from C x .^' by cyclic permutatins of X , Y , 
and Z factors in the right hand side of Eq. (|57j). In these notations, 

^ = E[ c cc a c^c + c cc^c«J] 
c 
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+ ^(C Cf ,a ( + C+,a+)(2n c , + 1) (58) 



with C^*, 



C^, — {Uj X r iUjy>)Cj.££l , 

C CC = ( U 3*' - iU iv') C r,tC' ■ ( 59 ) 

To calculate the probability of the spin-flip transition caused by a perpendicular-to-plane 
field E (t) || z , we express the z -coordinate operator in terms of operators (a^, a+) 

* = E k(k + ^), (60) 

where are defined in Eq. (jHSJ)- Substituting (JH%|) and (jHOj) into Eqs. (j%2} and (jSHJ), one 

arrives at 

L z {n^n v ) = (n^,^,^ \z BO \n^,n n ,l) 

= 4E(2«C' + 1)^- (61) 

with z/ = ± . Remarkably, despite the fact that diagonal and nondiagonal coefficients C^> 
enter into /C in a nonsymmetrical way, in (j61|) the symmetry is restored, and £ and £' take 
both values, (,(' = £,iV ■ Equation (|6ip can be conveniently rewritten as 

L z (n^n v ) = -7^k( 2,l (' + 1 ) 

x L%±% + . s %^%). (62) 

It is important for following calculations that the dependence of coefficients Cj.^i on 
indices ( and (' can be factorized (see Appendix |Bj) as 

C j M=p c ,[d j R K + {R < x D)j (63) 

with 

(3^ = —mu^/2h, P n = mu> rj /2h. (64) 

Therefore, the dependence of the coefficients Cj-^r on the last subscript is universal and 
rather trivial. Their dependence on the two first subscripts and the angles (6, 0) can be 
written in terms of three 3D vectors 

R ( = (X ( ,Y C ,Z C ), d=(b Y _ z ,b z ^ x ,b x „ Y ), 

D = (byz,bzx,bxY), (65) 
27 



of which the first one is defined by Eqs. (|21jh while d and D depend only on (8,<f)) . 
This structure of coefficients Cj-^c can be established by direct calculations. Formulas for 
coefficients b with different indices are given by Eq. (jB2|) . 

Details related to calculating matrix elements of Eq. (|62j) are described in Appendix [BJ 
We present the final expression following from Eqs. (jBf f|) in the form that is similar to 
Eq. ()44|) to simplify the comparison with the results of Sec. 1111 Al 

af(n i ,n ri )uj s u c sm9 . 
L z (nt,nr,) = ^ ^^-sign{a;o-w c } 

x [sin 2(f) cos 9 F s — i cos 20 F c ] , (66) 

where 

af{n^n n ) = a* D cu [to v (l + 2n v ) - c^(l + 2n ( )]/tt 2 (9) 

with a J, = —5muJo/2h. (67) 

Functions F s and F c are defined as 

,2cos 2 #-sin 2 



F s (8) = {u c - lu s 







2u c ujg — 3u>l 



1 + 3 sin 2 9 



(68) 



F c (6) = (u c cos 2 9 - u s ) + \{uo c uo 2 s (l + cos 2 9) 

-ul cos 2 9 + cj 2 ^ cos 29] . (69) 

It is seen from Eq. (JBTj) that the effective coupling constant of§ depends on the population 
of the both £ and rj - type energy levels. This factor reduces to a constant when (i) only the 
lower level of the spatial quantization is populated, n n = , and (ii) the Fermi energy is small 
compared to the spatial quantization energy, n^uo^ <C io n . In the strong confinement limit, 
with Q 2 w cg>q , we get aff ~ —5muo/2h . Having in mind that mu)o/2h w {^/d) 2 , d playing 
the role of the confinement length, we arrive at the usual estimate for c*d . However, when 
the criterion of strong spatial quantization is not fulfilled, cuff influences the dependence of 
the EDSR intensity on the polar angle of B . 

Functions F s and F c are defined in such a way that their first terms dominate in the 
strong confinement limit, ujq ^> uj c ,uj s , and these coefficients coincide exactly with the 
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corresponding coefficients in Eq. ()44|) . Therefore, in this case Eq. ()66|) reduces to Eq. ()44|) . 
However, it will be shown in the next section that when ujq ~ lu c , u s , the factors F s and 
F c influence profoundly the intensity of EDSR. 

A. Angular dependence of EDSR 

In Fig.EK,b the dependence of the EDSR intensity on the polar angle 9 has been presented 
for the Hamiltonian Hr describing SIA. It has been discussed at the end of Sec. IIII Al that 
the applicability of the Hamiltonian Hp is more restricted as compared to Hr . Therefore, 
for calculating the angular dependence of EDSR intensity coming from BIA in a wide range 
of polar angles 9 , we apply the general results derived in the previous section. Comparison 
of Eqs. (}4"4"j) and shows that this generalization influences EDSR intensity in three ways. 
First, through new terms in functions F s {9) and F c (9) . Second, through 9 dependence 
of Q(9) . Third, the transition probability acquires a dependence on population numbers 
n^, n v entering into the effective coupling constant a e ^(n^, n v ) . Also, as distinct from Hr , 
the Hamiltonian 7i,D does not possess continuous rotational symmetry about the z axis. 
Therefore, we calculate the 9 -dependence of the EDSR intensity for two specific values of 
0, = and = 7r/4. 

The polar angle dependences of the EDSR intensities presented in Fig. |H1 have been 
calculated for the same values of uj c /ujq and lu s /uq as the data of Fig. Eland n e = 2ul(9 = 
0) . We have also used the same Gaussian level broadening. Because we disregarded the effect 
of SO coupling on the energy spectrum, all the above discussion of the energy spectrum and 
the 9 -dependence of the chemical potential rj{9) is completely applicable to the current 
case. Hence, the basic behavior of the intensity including the principal maximum near 
co>£ = \u s \ and existence of a system of up- and down-cusps remains intact. The difference 
comes from the envelope factors F s (9) and F c (9) , from the # -dependence of a e § , and from 
the fact that for BIA the matrix element of EDSR does not vanish when 9 — > 7r/2 . For large 
9 the cyclotron frequency uo c {9) = u c cos9 vanishes and our theory becomes inapplicable. 
Therefore, we restrict ourselves with the region 9 < 0.9 x n/2 . 

In Fig. EH, the EDSR intensity is displayed for = 0. It shows cusps similar to Fig. Ek, 
however, all features seem small in the scale of the figure. This difference in the scale can 
be understood from the shape of the envelope function shown in the inset. It includes the 
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factors from Eq. (|66|) that depend on 9 explicitly and do not involve quantum numbers 
n^,n v . As distinct from the envelope function of Fig. that includes a factor sin 2 29 , 
Eq. (|43|). this envelope function changes slowly on the right from the principal maximum at 
td£ ~ \ou s \ and does not vanish at 9 = n/2. Therefore, the EDSR intensity remains strong 
for large 9 , and we have shown a curve for a very small V = 0.02/^0 to demonstrate the 
persistence of the principal maximum. At larger T values the maximum washes out, and 
for uj c cos9 < T/h our approach is no more applicable. Remarkably, for u c = Uq/2 the 
envelope function is close to its shape in the strong confinement limit as one can see from 
the inset. 

In Fig. \§jp, the intensity of EDSR is shown for the same value uj c /u = 0.5 but for 
<p = 7r/4. In this case the envelope function includes a factor sin 2 29 and vanishes for 
9 = tt/2 . For this reason, Fig. bears much more similarity with Fig. Deviation of the 
envelope function from its strong confinement limit is by a factor of about 2. Comparison of 
Figs. and Hb suggests a considerable azimuth dependence of EDSR intensity, especially 
in the large 9 region. 

It is very instructive to compare panels c and d of Fig. drawn for u c / uoq = 2 with 
the panels a and b of the same figure. The first difference is easily seen from the envelope 
functions shown in inserts. They show zeros for intermediate values of 9 , 9 « 0.82 x 7r/2 
and 9 « 0.61 x n/2 for = and = vr/4, respectively. The corresponding minima 
originate from the zeros of functions F c (9) and F s (9) , respectively, appear only for specific 
values of the azimuth , and strongly influence the EDSR intensity. 6a The existence of an 
additional zero that is common for both panels c and d (but absent from a and b) and 
persists for arbitrary is the second difference. It originates from the zero of a^in^^n^) 
at 3uj v (9) = oj^{9) that is located at 9 w 0.42 x n/2 . For this 9 value, the filling factor v 
is within 2 < v < 3 . Hence, the spin flip transitions involve the = 0,n v = 1 electrons. 
The third important difference is a strong enhancement of EDSR intensity at u c = 2ujq as 
compared with u c — Uq/2 . Intensities in Fig. IHlare given in arbitrary units, but these units 
are the same for all panels when normalized on the confinement frequency ujq . The origin 
of the enhancement can be easily understood from Eq. (|4^|) that in the strong confinement 
regime includes a factor ar,uo s /ujQ with oc Uq . Therefore, when B increases, one can 
expect enhancement in EDSR intensity by a factor of about (uj s /ujq) 2 . 
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V. DISCUSSION 



The results of our calculations allow one to evaluate the intensity of EDSR and to compare 
it with the intensity of the usual paramagnetic resonance (EPR) caused by a time dependent 
magnetic field. Their ratio depends on the SO coupling mechanism and the geometry of the 
experiment. In both the EDSR and EPR experiments, the efficiency is controlled by two 
major factors: (i) by the characteristic length I which in the case of EDSR is equal to the 
matrix element of the spin-flip transition of Eq. (|33J). and (ii) by the population difference 
of two spin sublevels that depends on the 2D electron concentration and the polar angle 9 
between the external magnetic field and direction of 2D confinement. 

The effect of the latter factor does not depend on the mechanism of spin-flip transitions 
and influences in a similar way both the EDSR and EPR. It manifests itself in the modula- 
tion of the spin-flip transition intensity as a function of 9 and depends on the population of 
spin sublevels controlled by the 9 dependent position of the chemical potential rj(9) . The 
corresponding up- and down-cusps are seen in Figs. El and El showing the polar angle depen- 
dence of EDSR intensity for both basic SO coupling mechanisms. Level broadening washes 
out this population dependent fine structure. However, for a reasonably small broadening, 
r/hu c ^ 0.1 , this fine structure should allow studying the energy spectra of quantum wells. 

A geometry with the electric field E perpendicular to the confinement plane, similar to 
the one used by Kato et alM, in principle allows access to electron spins at a nanometer scale. 
However, spin-flip transitions occur only in a tilted magnetic field and vanish when 9 = 0, 
as is seen from Eqs. (JH3), (HU), and (JHHJ). Indeed, it is a tilted magnetic field that couples 
the in-plane and perpendicular to the plane electron motions and allows a perpendicular to 
the plane field E to produce cyclotron and spin-flip transitions. 

The efficiency of the Dresselhaus mechanism of EDSR in this geometry can be evaluated 
by using the effective length = L z {n^ = 0,12^ = 0) defined by Eq. In the strong 

confinement limit, c^o 3> uj c ,uj s , this length is about ~ (a* D /huo)(uJ s /Ld c )(u) c /uJo) when 
expressed in terms of the effective coupling constant a* D = —8muJo/2h of Eq. (Jfi7Jl . The 
factor uj c /ujq reflects the fact that the deviation of the system from a strictly 2D geometry, 
ujq — > oo , is critical for the gate- voltage controlled EDSR. The numerical factor uj s /uo c = 
gm/2m is about 0.16 for GaSb and InAs and about 0.34 for InSb. Even for a weak 
magnetic field oj c /ujq = 0.1, we estimate ~ 5 x 10 -11 cm to 5 x 10~ 10 cm using a 
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typical value m ~ 0.05mo for the mass and also 5 ~ 20 eV A 3 for GaAs and 200 eV A 3 
for InSb and GaSb. 64 Hence, under these least favorable conditions the EDSR length is 
comparable to the similar length for EPR that is about Zepr ~ \g\ ~ 10~ 10 cm, with 

— h/niQC Ri4x 10 -11 cm for the electron Compton length and \g\ ~ 10 for the electron 
g factor. 

Remarkably, there are two factors related to the electronic confinement in a quantum well 
and the experimental geometry that can increase Ijj essentially. The first one is related to 
the resonance behavior of lp in the angle range near u c cos8 ~ \u s \ . The resonance in 1^(9) 
is cut off by the level broadening, and the increase in the effective length is about a large 
factor of uj c /T ^> 1 . The second factor is related to the ratio u> c /ujo reflecting the deviation 
of the system from the strict 2D geometry; it is not necessary small. E.g., in the Kato 
et al. experiment^ this ratio was about u c /uq ~ 0.5 because a wide parabolic well with 
an effective width about 50 nm and a strong magnetic field B = 6T were used. A nearly 
sixteen time increase of the EDSR intensity due to a four time increase of uj c /ujq from 0.5 
to 2 is seen in Fig. where the intensity of EDSR for all cases was calculated in arbitrary 
but the same units based on a fixed confinement energy ujq . The increase in uo c /uoq leads 
to a similar increase in the intensity of the EDSR controlled by Rashba spin-orbit coupling 
as one can see from Fig. El 

The magnitude of the Dresselhaus spin orbit coupling constant a* D is controlled by 
the bulk parameter 6 and cannot be modified significantly. Even in the wide quantum 
wells used in the Kato et al. experiments,— the effective 2D coupling constant reaches 
only the value a* D ~ 0.3 x 10~ 10 eVcm. This is much less than the typical value of the 
Rashba constant an ~ 10 -9 eVcm for InAs based quantum wells.— Even larger values an m 
(3-=-6) x 10 -9 eVcm were reported in Refs. 23j and jf^]. Large values should significantly 
enhance the efficiency of EDSR whose magnitude can be estimated using Eq. P3*|) and gives 
the effective magnetic length lj^ ~ (an/tiuJo)(uj s /uj c )(uj c /uo) . This suggests that asymmetric 
quantum wells should have a potential for enhancing EDSR. However, calculations of an 
depend strongly on the boundary conditions,-^ and the dependence of cxr on the well width 
has not been investigated yet. An increase in ujq might happen to be a price for increasing 
an . From this standpoint, the potential of wurtzite type materials where an emerges as a 
bulk parameter should be investigated in more detail. 

The estimates provided above show that electrically manipulating electron spins is prefer- 
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able to the magnetic one not only because it allows access to electron spins at a nanometer 
scale but also because a significantly larger coupling constant can be achieved. In order to 
reach high efficiency of EDSR in the geometry with an electric field E perpendicular to 
the quantum well plane, one needs to use wide asymmetric wells and a strong magnetic field 
B strongly tilted to the quantum well plane. The EDSR should be also stronger in narrow 
gap semiconductors where the factor uj s /uo c is typically larger. 

Let us estimate now the advantages of EDSR with an in-plane electric field E . The ef- 
fective lengths characteristic for the Rashba and Dresselhaus spin-orbit coupling mechanisms 
are described by Eqs. flHJ) and (JH2J), respectively. By the order of magnitude, V R ~ a R /hu c 
and 4 ~ ai£)/huj c for B || z . A comparison with the estimates used above for the ge- 
ometry with a perpendicular to the plane field E shows that both small factors in the 
expressions for the effective lengths do not appear in the equations for the in-plane geome- 
try. The small factor uj c /uj has been discussed above and is related to the deviation from 
the strict 2D limit. This criterion, as well as necessity of using a tilted magnetic field B are 
specific for the geometry with a perpendicular to the plane E and are required to couple 
this field to in-plane dynamics. With an in-plane electric field the both limitations do not 
exist any more, and this critical difference manifests itself in Eqs. ()49|) and (|52|) . 

The second small factor, uj s /ujo, originated because of the confinement in the electric 
field direction; a similar factor appears for EDSR with impurity centers 6 . For an in-plane 
electric field, the motion in the direction of E is unrestricted and that is why the parameter 
uj s /uJo does not appear in Eqs. (jHH) and (JH2*j) . 

With typical values of old ~ 0.3 x 1CT 10 eVcm and a R ~ 10 -9 eV cm, m ~ 0.05mo, 
and B « 1 T, we get l D ~ 0.3 x 1CT 6 cm and l R ~ 1CT 5 cm correspondingly. As a result, 
4> 4 3> %, l R as well as if,, l R 3> Zepr , and electrical spin operation by an in-plane electric 
field should be especially efficient. Recent experimental data^ are in agreement with this 
prediction of our theory.— 

The most important quantity characterizing spin operation efficiency by a resonant elec- 
tric field E(t) is the Rabi frequency Vt R = eEl/h. With / w 4,4 ~ 10~ 5 -r- 10~ 6 cm 
as estimated above, we find that fl R ~ 10 10 s _1 in electric fields as small as only about 
E ta 0.6 4- 6 V/cm. 

This extraordinary efficiency of in-plane electric fields bears a promise of operating spins 
by a set of small vertical gates producing considerable in-plane fields between them. Such 
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a geometry can allow combining the access to electron spins at a nanometer scale typical of 
vertical gates with the high EDSR intensity typical of in-plane electric fields. 

In conclusion, our theory predicts efficient EDSR with free electrons in quantum wells. 
Usually, its intensity should be much higher than the intensity of the standard EPR. Elec- 
trical spin operation should be extraordinary efficient with an in-plane electric field. In the 
perpendicular to the plane electric field geometry, EDSR requires mixing of the in-plane and 
perpendicular to the plane orbital motions which is provided by a tilted magnetic field. As 
a result, EDSR in this geometry is significantly weaker than in the in-plane geometry, and 
semiconductor compounds with large g -factors are highly advantageous. The dependence 
of the EDSR intensity on the magnetic field direction, on both its polar angle and azimuth, 
is indicative of the competing mechanisms of spin-orbit coupling. 
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APPENDIX A: USEFUL IDENTITIES 

All formulae for intensities of EDSR, as well as for eigenfrequencies uo^ and uu v , include 
the auxiliary angle 7 that should be found from Eq. ©• However, it can be usually 
eliminated and equations for the intensities of EDSR can be simplified due to the existence 
of a system of identities that will be derived in this Appendix. 

It follows directly from Eqs. ()10|) that 



2,2 2,2 
U)£+U v = uj c +uj . 



(Al) 



Using Eqs. and (fTUj). one arrives at the identities 



22 22 2 n 

UOcUOr, = UJ r .UO n cos U 



(A2) 



cu 2 sin 2 7 + uJi cos 2 7 = u)q cos 




(A3) 
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£j| cos 2 (fl + 7) + sin 2 (fl + 7) = < cos' 0. (A4) 
^From Eqs. (|A2|) and ()A3|) follows the identity 

sin 2 7 cos 2 7 1 

— + — ^ = — A5 

u 2 v u 2 u 2 c 

that ensures fulfillment of the commutation relations [k x , x\ = [k y , y] = —i . Subtracting 
Eqs. (HH and J5H) from (EH), one finds 

22,2-2 2,2-2/1 

u>g cos 7 + a> sm 7 = uj c + uj sin 0, (Ao) 
w| sin 2 (# + 7) + w 2 cos 2 (fl + 7) = uj 2 sin 2 + wjj. (A7) 



2 



Using Eq. ()A7|) . one immediately finds 



cos 2 (fl + 7 ) sin 2 (fl + 7 ) ^ 2 + ^ 2 sin 2 fl-^ 

2 ,2 + „,2_,,,2 („,2 . .3 , .7\ ' 



^|-^ S 2 )K-^ S 2 ) 



and applying Eqs. ()A2|) and ()A4|) . one arrives at 



cos 2 (fl + 7 H 2 sin 2 (g + 7 V 2 _ o; e 2 cos 2 ^K 2 -^ 2 ) 

Two similar equations can be written for the angle 7 . Remarkably, the right hand sides 
of Eqs. (jA3|) - (jA7|) do not depend on 7 , and the same is true for the right hand sides of 
Eqs. flA"8j) and (lA"9|) because 

(a; 2 - u; 2 )(a; 2 - u 2 s ) = ulu 2 c cos 2 - oo 2 s (uj 2 + uo\ - oo 2 s ) } (A10) 

as follows from Eqs. (jAl|) and ()A2|) . 

^From Eqs. © and ()10|) one can also derive an identity 



00 



|tan 7 = w 2 tan(^ + 7) , (All) 



from which the following relations follow 

sin(6 l + 7) cos 7 shi7cos(# + 7) sin 27 



oo 2 oo 2 loo 2 cos 6 



sin 2 (# + 7) 
2c<j2 cos 6* 



(A12) 
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In turn, a set of new identities follows from Eq. (|A12|) 



co>| cos(# + 7) ^COS7 



COS7 = -| — , cos(6' + 7) 



UJ 2 COS 9 UJq cos 9 



u 2 sin(fl + 7) . w|sin7 

sin 7 = - — , sm(^ + 7) = -| -. 

uJc cos 9 Uq cos 9 

(A13) 

These identities, together with Eq. (|2Tj). allow proving the commutation relation [x, y] = 0. 
Also, using Eq. (JHJ), one easily finds 

sin 29 

tan2 7 = - — — — , 

(uj c /uj ) z — cos 29 

, . sin 29 . 

tan 2 + 7 = — — — — . A14 

cos 20 — (uj /Lj c y 

These equations allow one to derive the identity 

o n Unsm29 ur?sin 26* , , 

- ^ = = — ^— . (A15) 

c " sin 2 7 sin2(# + 7) v ; 

All the above identities have been derived directly from Eqs. © and (fTU|) and are valid 
for an arbitrary ratio uo c / u and an arbitrary value of 9 , < 9 < 71 . 

Finally, using Eqs. (jA14|) . one can derive explicit expressions for sin 27 and sin 2(9 + 7) 
that allow one to find and oj v from Eqs. (|A1|) and (|A15|) . For < 9 < %/2 , 

sin 27 = -— I sin26 l sign{^ - co c } , 
u 2 

sin 2(9 + 7) = sin 20 sign{cu - w c } , (A16) 

with f2 2 (#) defined by Eq. (j!5|) . Equation ()A14|) defines 7 with the accuracy to the phase 
±7r/2 . The choice of the phase made in Eq. (jA16|) results in the energy spectrum of Eq. ()14|) 
shown in Fig. ^ With this choice of the phase, Eq. (|A15j) takes the form 



OJ 



2 v -Lul = n 2 (9)sign{LU -uj c }. (A17) 



It indicates that the difference cjg — uj v , when considered as a function of u c / ujq , changes 
sign at u c = ujq . 



APPENDIX B: SEVERAL SUMS CONTRIBUTING TO EDSR 



For the calculation of Cj.^t we need to find bilinear combinations of the coefficients 
Xg , Y^t , and appearing in the brackets of Eq. |37| Straightforward calculation of these 
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coefficients based on their definition by Eq. (|21|) results in 

|X ? | 2 -|r c | 2 = Pcb X -Y,X ( Y£ + Y c X2 = fcb XY , 

and b YX = b XY , (Bl) 

where (3^ = -mw^/2/i, j3 v = moo v /2h, and all other combinations can be found from 
Eq. (jBl|) by cyclic transpositions of X , Y , and Z . Here 

bx-Y = - sin 2(9 + 7) tan 9 cos 20 , 

b Y -z = -[cos2(6» + 7) + isin2(6' + 7)tan6'cos 2 0] , 

b z _ x = [cos2(6» + 7) + i sin 2(6* + 7) tan 6^ sin 2 0] , 

bxY = - sin 2(8 + 7) tan 9 sin 20 , 
byz = sin 2(9 + 7) sin , 

b zx = sin 2(fl + 7 ) cos . (B2) 
Sums standing in the numerators of Eq. (jfilj) involve the expressions 



Rj£ + R*£ = 2cu c y m/2hui£ cos 7 G\- 



Rjr, + R* jv = 2cu cV /m/2^sin7 GJ, (B3) 
and differences standing in the same numerators involve the expressions 









p* 




p* 


= (- 


sin 0, 1 



2iJmu ( /2h f ( G~ 



-2iJmu c /2h f c G~ 



(B4) 

where 

G + = (— sin0, cos0, 0), G ~ = (cos0, sin 0, 0), (B5) 

C is defined as £ = 77 and 77 = £ , and is defined as = 0, m,, = 1 . 

Using Eqs. (jB3|) and (jB4|) . it is convenient to rewrite two groups of coefficients needed 
for calculating numerators of Eqs. (JoTj) as 

/ s ^sin2(0 + 7) . , 
l&dCjM + Ck,<>) = -(-ir " 2uJ CQs6 - Pt'Mt, (B6) 
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+ (1- 



i(-irc sin 2(^ + 7 )M; 



(B7) 



where M and A?" are 3D vectors 



Mf = Gfdj + (G ± x D 



N = (b Z x, -b Z Y, -bx-r) ■ 



(B8) 



Using these expressions, one can transform Eq. (jfilj) by collecting terms having common 
factors 2n^ + 1 . Then the following sums emerge 



C rXC + C j 



zC^C ,.,2 2 



+ /? c a; c a;;?sinfl 
j V(6) 



(B9) 



< , ,2 _ , ,2 



'2 J ' Z>(0) j X>(0) 



Plugging these equations into Eq. (fB2~jl. we arrive at 



££(2n f + l)/? c 



Z>(0) 



1 , . cUoU; 2 sin 26 1 , .cjJu; 2 cos 2 # — cj 2 ) 
— — — — h iM — 

2 D(0) + D(0) 



with 



^± = ( r. v - iU jy ,)Mf, AT = (C,y - iU ry> )N j 



Calculating these coefficients by using Eqs. (JHJ) and (jB2|) . we find 



(BIO) 



(Bll) 



sin 20[(7/4) sin 2(5 + 7) sin 5 + cos 2(5 + 7) cos 5] 
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— % cos 20 cos 2 {9 + 7) , 
M- 

3 

= - cos 20[cos 6 1 cos 2(6> + 7) + - sin 2(6* + 7) sin 9] 

+ i sin 20[^ sin 2(0 + 7) tan 9 - cos 2(0 + 7)] , 
J\f 

= sin 2 (0 + 7) [cos 20(cos + — sin 9 tan 9) + i sin 20] . 

(B12) 

From these equations follows the final equation (jfifij) . 
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FIG. 1: (Color online) The energy spectrum of electrons confined in a parabolic quantum well 
in a tilted magnetic field. Here, ujq is the confinement frequency, tu c = eB/mc is the cyclotron 
frequency, 9 is the polar angle of the field B , uj^ and uj v are the frequencies of two eigenmodes 
, C = £>?7 ■ ( a ) The dependence of uj^ on the ratio oo c /^o f° r three values of 9 . Cyclotron-like 
modes are shown by dotted lines and confinement-like modes by full lines. For 9 = these modes 
intersect, while for 9^0 they interchange at u c = uq . (b) The dependence of oo^ on 9 for 
three values of uj c /ujo . The modes retain their identity in the whole interval < 9 < ir/2 , and 
(jj T] > uj^ for lu c < ujq while uj T] < uj^ for oj c > ujq . At the degeneracy point, uj c = ojq , the two 
branches have no specific identity. They are separatrices dividing the plane into one rj (gray) 
and two £ (white) regions. 

FIG. 2: (Color online) The difference of the filling factors, Ais , of the two components of a spin 
doublet of an isolated energy level as a function of the total filling factor of this level < v < 2 . 
For a small level width T , its shape is close to a triangular form, and an up-cusp develops at 
v = 1 . In a plot including adjacent energy levels, down-cusps at the end of this interval will arise. 
With increasing T, the cusps are smeared. However, at T as large as T = huj s /2, pronounced 
extrema at v = 1 and v = 0,2 are still seen. 
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FIG. 3: (Color online) The dependence of the EDSR intensity on the polar angle 9 of a magnetic 
field B for Rashba spin-orbit coupling, Eq. (j3J). The electric field is perpendicular to the confine- 
ment plane, E || z . The ratio of spin-flip and cyclotron frequencies, lo s /lo c = —0.17, is typical 
of InAs; ujq being the confinement frequency. The calculations were performed for three values of 
the level width T , as specified in the figure. For B || z , electrons completely populate the ground 
level, = n v = . (a) Strong confinement, lo c /u>q = 0.5 . With increasing 6 electrons populate 
higher cyclotron-like £ levels, while they remain at the ground r] level, (b) Weak confinement, 
u) c /u)q = 2 . With increasing 6 electrons populate higher r] levels, while they remain at the ground 
£ level. In both cases, the principal peak originates from the resonance uj^ = \uj s \ with C = £ f° r 
(a) and C = V f° r (b). The fine structure of the spectrum, distinctly seen for small level width 
r = 0.05ujq , originates from populating higher levels in a tilted magnetic field B . Up-cusps 
are marked by up-arrows and down-cusps by down-arrows; for details see text. In (b), the inflexion 
point seen between the v = 5 cusp and the EDSR maximum reflects an abrupt change in the pop- 
ulation of intersecting levels at oo^ = \uj s \ . Inserts show the dependence of the chemical potential 
r]{v) on the filling factor u and the quantum numbers of partially populated states, or , 
that contribute to spin-flip transitions inside the various regions of filling factor values. 



FIG. 4: The dependence of the EDSR intensity on the magnetic field direction for a (001) quantum 
well. The electric field is perpendicular to the QW plane, E || z . The electron concentration is 
low , hence, for all 9 values electrons populate only the ground level, ng, = . The level width, 
r = 0.2a;o ; oj c ju)Q = 0.5, uj s /u c = —0.17. (a) 2D Dresselhaus spin-orbit coupling mechanism, 
Eq. |2J). The four- fold symmetry of the figure reflects the symmetry of the square of the spin- 
flip transition matrix element, (b) Interference of 2D Dresselhaus and Rashba spin-orbit coupling 
mechanisms. The two-fold symmetry of the figure reflects the C 2v symmetry of the Hamiltonian 
Hp ; coupling constants have been chosen as old = olr . 
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FIG. 5: The dependence of the EDSR intensity on the magnetic field direction for a (001) quantum 
well. The electric field is in the QW plane, E _L z . The electron concentration is low; parameter 
values are the same as in Fig. (a) and (b) - Rashba spin-orbit coupling, Eq. @. For (a,) the 
electric field E is parallel to the projection B j_ of the field B onto the confinement plane, 
E || B x , hence, the intensity is rotationally symmetric. For (b), the electric field is parallel to 
the (110) crystallographic axis, E || (110) . (c) and (d) - 2D Dresselhaus spin-orbit coupling, 
Eq. (j3J). For (c), the electric field E is parallel to B ± , E \\ B ± , hence, the intensity shows 
four- fold symmetry. For (d), E || (110) , and the gross features are similar to those of figure 
(b). For (e), E \\ B±, the two spin-orbit coupling mechanisms with qd = ctR interfere; the 
symmetry is lowered because of their interference. 



FIG. 6: (Color online) The dependence of the EDSR intensity on the polar angle 9 of a magnetic 
field B for Dresselhaus spin-orbit coupling, Eq. @ , for two values of the azimuth <j) . The electric 
field is perpendicular to the QW plane, E || z . The notations and basic parameter values are the 
same as in Fig. |3 Similarly to Fig. |3J the spectra show the principal maximum at = \u s \ and 
the fine structure of up- and down-cusps indicating repopulation of different energy levels with 9 
changing. The major distinctions from Fig. include (i) a minimum at 9 = 0.42 x tt/2 originating 
from a zero in the effective coupling constant a c ^ (n^,n v ) of Eq. 1)6 7|) . (ii) additional minima in 
panels c and d originating from the zeros of functions F c {6) and F s (9) at 9 = 0.82 x tt/2 and 
9 = 0.61 x 7r/2, respectively, and (iii) smearing of the main peak in panel (a) because of the 
slow angular dependence of the envelope function shown in the inset to that panel. For details 
see the text. The inserts to panels (a) and (b) demonstrate that the results found for the strong 
confinement limit still retain reasonable accuracy for uj c = ujq/2 . The inserts to panels (c) and (d) 
demonstrate the positions of the zeros of functions F c and F s . 
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